Early Higgs Boson Discovery in Non-minimal Higgs Sectors 



O 
(N 



Oh 
I 

Oh 



> 

00 
On 

m 

O 



X 



Spencer Chang,^^^'^ Jared A. Evans/'l^and Markus A. Luty^'l^ 

'Physics Department, University of California Davis 

Davis, California 95616 

'^Department of Physics, University of Oregon, Eugene, OR 974-03 

Particle physics models with more than one Higgs boson occur in many frameworks for physics 
beyond the standard model, including supersymmetry, technicolor, composite Higgs, and "little 
Higgs" models. If the Higgs sector contains couplings stronger than electroweak gauge couplings, 
there will be heavy Higgs particles that decay to lighter Higgs particles plus heavy particles such as 
W, Z, and t. This motivates searches for final states involving multiple W, Z, t, and bb pairs. A two 
Higgs doublet model with custodial symmetry is a useful simplified model to describe many of these 
signals. The model can be parameterized by the physical Higgs masses and the mixing angles a and 
P, so discovery or exclusion in this parameter space has a straightforward physical interpretation. 
We illustrate this with a detailed analysis of the process gg ^>- A followed by yl — > hZ and h — > WW. 
For niA — 330 GeV, mh — 200 GeV we can get a 4.5ff signal with 1 fb~^ of integrated luminosity 
at the Large Hadron Collider. 

PACS numbers: 



Introduction — The Large Hadron Collider (LHC) is 
currently exploring the high-energy frontier at the TeV 
scale, and its main goal is to discover the origin of elec- 
troweak symmetry breaking. Most of the effort in this 
direction is devoted to the search for a Higgs boson with 
mass below roughly 150 GeV. This is motivated by the 
fact that the standard model with a light Higgs boson 
gives a good fit to precision electroweak data, and by su- 
persymmetric models that predict a light Higgs boson. 
However, the minimal supersymmetric model generically 
predicts m^ < mz = 91 GeV for the lighest Higgs mass, 
violating the experimental limit m^ > 114 GeV from 
searches at LEP [T]. The Higgs mass can get loop con- 
tributions that avoid the experimental bound, but at the 
cost of fine-tuning UV contributions to the Higgs mass at 
the percent level. Eliminating precisely this kind of tun- 
ing is the primary motivation for supersymmetry, and 
this has motivated a great deal of work on supersym- 
metric models to reduce this tuning (see for e.g. [SHI] 
and [5^ for a more complete list of references). With or 
without supersymmetry, the precision electroweak fit is 
compatible with a heavier Higgs boson if there are ad- 
ditional particles with masses and couplings that break 
electroweak symmetry (see for e.g. [Sj). 

In this Letter, we study the phenomenology of non- 
minimal Higgs sectors with Higgs boson masses above 
2mw. Non-minimal Higgs sectors occur in many frame- 
works for avoiding fine tuning of the Higgs sector, in- 
cluding new strong dynamics ("technicolor") (for a re- 
view, see f7|), composite Higgs (see for e.g. [Hlin]), and 
"little Higgs" theories (for reviews, see [lUllII]). Heav- 
ier Higgs bosons are associated with Higgs sectors with 
stronger self-couplings, and are therefore particularly 
well-motivated in technicolor and composite Higgs mod- 
els. In particular, Ref. [12] argued that spin-0 resonances 
are a generic and prominent feature of models with strong 



dynamics in the electroweak symmetry breaking sector, 
and the present work originated in the construction of 
simplified models to describe their collider phenomenol- 
ogy [13] . However, the collider signals we discuss are 
applicable to a much wider range of models. 

Higgs bosons couplings are generally proportional to 
the masses of the particles to which they couple, so the 
Higgs particles tend to decay to the heaviest particles 
that are kinematically allowed. In the present context, 
this implies decays to lighter Higgs bosons, t, Z, W and 
b particles. In particular, we can have cascade decays 
leading to states containing several heavy standard model 
particles. 

Simplified Model — To explore the phenomenology, we 
use a 2-Higgs doublet model as a simpHfied model [l4] . 
This model possesses built-in, good high-energy behav- 
ior, has a simple parameter space, and can be simulated 
using state-of-the art Monte Carlo tools. We impose 
SU{2) custodial symmetry as well as a discrete symme- 
try H2 i-T- —H2 on the potential to simplify the model. 
The Higgs potential is then 



V 



mlH\Hi +mlHlH2 



+X^{HlHi){HlH2 



'-{HIH2 
A4 



(1) 



{HIH2 + h.c.)^ 



The physical particles consist of the CP-even neutral 
scalars ft,, H and a degenerate custodial SU{2) triplet 
[A.H^). The 6 parameters in the Higgs sector can be 
taken to be w = -y/uj -I- w| = 246 GeV, ttt,^, mf{, rriA — 
mH± , and the angles a and (3 defined conventionally by 
tan/3 = ui/wa, h = [cosaRe{H^) - suiaRe{H^)]/V2. 
The custodial SU{2) symmetry allows us to choose Hi 
to be the field that couples to the top quark, so all phe- 
nomenologically important couplings of the Higgs bosons 
are determined by these parameters. 



Phenomenology — We now discuss the signals beyond 
those of the standard model when some of the Higgs 
bosons are heavy. If the A is relatively light and has 
unsuppressed coupling to the top quark, it can be pro- 
duced via gg ^ A (via a top loop) followed hy A ^ tt 
or hZ. If the A is light but has suppressed coupling to 
the top quark (large tan /?) we have gg ^ H followed by 
H — >• AZ — ^ hZZ. If the A is very heavy, the domi- 
nant signal beyond the standard model is gg -^ H ^^ hh. 
All of these signals produce light Higgses, h, which de- 
cay either to bb (for rrih < 2mw) or WW/ZZ (for larger 
mfi). There are many modes, with the common feature 
that they involve production of multiple heavy standard 
model particles. 

In this work we study the process gg — ^ A^*^ — > Zh. 
If mh < 2mw the dominant decay is /i — > bb; in the 
context of supersymmetry, this has been discussed as a 
discovery mode of two Higgs bosons [15] . Zh production 
also occurs in the standard model, and backgrounds can 
be suppressed by focusing on the kinematic region where 
the h is boosted [111 HH- In the present scenario this 
is enhanced by the intermediate A in the production. 
Requiring prih) > 200 GeV the intermediate A is off- 
shell in our model, but there is still an enhancement of 
10 times or more compared to the standard model, as 
illustrated in Fig. [T] Since this signal is already under 
detailed investigation by the LHC experiments, we will 
not discuss it further here. 

We will study the case of a heavier h, where we have 
gg — > Zh — > W^W^Z via an on-shell A. To our knowl- 
edge, this has not been previously investigated in the 
literature. We focus on Z ^ l'^l~ (/ = e or ^), which 
suppresses all backgrounds that do not involve a Z. We 
consider WW — > jjl"^ + p.^, which has a 29% branching 
ratio and allows us to reconstruct both the h and the A. 
The final state WW — > l^l'^ + prp has branching ratio 
4.5% and is very clean, but does not allow reconstruction 
of the h (or A). Another possibility that we do not in- 
vestigate here is Zh — !• ZZZ -^ jjl^l^l^l^ , which has a 
branching ratio approximately 0.2 times the mode we are 
considering. With more data, these modes will be useful 
in confirming the signal and its interpretation. 

Because we require a leptonic Z, the only important 
backgrounds are those involving a Z. The largest back- 
ground is WZ +iets. Other backgrounds we consider are 
Z + jets, itZ, and ZZ + jets. The signal was generated 
with up to 1 additional jet using MadGraph [18] with the 
2HDM4TC model files [13], and applying a i^- factor of 2.4 
(extrapolating from [H]). Backgrounds were simulated 
using Alpgen [5D] with up to one additional jet. Jet show- 
ering was done by Pythia [2T] with MLM jet matching 
[22] for both the signal and backgrounds, and the de- 
tector response was simulated using PGS [23]. The K- 
factors used for the backgrounds are: Kyyz = 1.6 [Ml US], 
Kz = 1.6 [ig, Kzz = 1.3 [m HH], and Kuz = 1-3 
[25] . We study a benchmark model with niA = 330 GeV, 



e 

= 160 






^ 


^ 














^140 


:/' 






\ 














Sl20 

o 


r 






\ 














100 


- 






\ 














80 
60 
40 


r^ 






, 1 


V. 





:r 


1 , 


/ 


y'^" 

/ 


20 








300 



350 



400 



450 500 

A° Mass (GeV) 



FIG. 1: Cross-section times branching ratios at the 7 TeV 
LHC as a function of ttia. Solid line: pp — )■ A — )■ Zh — >■ 
W^W'Z -> l+ri^jj+^rj, for a model with ruh = 200 GeV, 
tan/3 = 1, sin a = 1. Long dashed line: pp — )■ A* — )■ Zh — >■ 
bb + Pt with pT{h) > 200 GeV and f:,^, > 200 GeV for a 
model with ruh = 120 GeV, tan/3 = 1, sin a = 0. The com- 
parable standard model cross section is 2.7 fb. Short dashed 
line: pp-^ A' -^ Zh^ l+rtb with pt(/i) > 200 GeV for the 
same model. The comparable standard model cross section is 
1.1 fb. 



irih = 200 GeV, mn — 1 TeV, sin a = 1, and tan/3 — 1. 
The cross section for producing A at the Tevatron is 
190 fb. The CDF experiment has searched for the WW Z 
final state from Zh associated production in the standard 
model |30j. Although this is a neural net analysis that 
cannot be directly compared with our model, the limit 
on the WWZ cross section is approximately 450 fb. We 
conclude that Tevatron searches are not senstive to this 
model, and we focus on the prospects for discovery at the 
7 TeV LHC. 

For an LHC search, we require the events to pass one 
of the following lepton triggers: (i) single lepton with 
prie) > 30 GeV or pt(m) > 20 GeV; (ii) double lepton 
with prih) > 17 GeV, pt(^2) > 10 GeV; (iii) triple lep- 
ton with pt{Ii,2.3) > 10 GeV. All triggered leptons are 
required to be central, \T]e\ < 2.4 and l??^! < 2.1. These 
triggers are similar to those used in the 2011 LHC run. 
The selection cuts then require at least 3 central leptons 
with pt(0 > 8 GeV and at least 2 jets with \gj\ < 2.5 
and ptU) > 30 GeV. We then impose the following 
cuts: (i) Px > 20 GeV. This effectively suppresses 
the Z + jets background, (m) We require that 2 same- 
fiavor, opposite-sign leptons reconstruct to the Z with 
Ifnii — mz\ < 7 GeV. We estimate that loosening this 
cut further will allow tt backgrounds to be significant. 
[iii) We combine missing px and the hardest remaining 
lepton (and, if that fails, the next hardest remaining lep- 
ton) to reconstruct to the W mass. This will give either 
2 solutions or no solutions, {iv) Two of the jets are re- 



Process 


Selection 


Pt 


Z{U) 


VKOPt) 


W(jj) 


WZ + 2,3 iets 


8.85 


7.94 


7.50 


5.70 


1.65 


tiZ + 0, 1 jets 


0.932 


0.896 


0.675 


0.475 


0.293 


^ + 2,3 jets 


2.30 


0.848 


0.587 


0.586 


0.226 


^2' + 1,2, 3 jets 


1.37 


0.611 


0.573 


0.504 


0.180 


Total Background 


13.4 


10.3 


9.34 


7.26 


2.35 


Signal +0, 1 jets 


29.8 


25.9 


22.3 


16.3 


9.89 



TABLE I: Signal and background cross sections (in fb) for the 
benchmark model at the 7 TeV LHC. The cuts are described 
in the text. The incorrect solutions from reconstructing the 
leptonic W are not counted in this table. 
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FIG. 2: Reconstructed h mass for 1 fb~^ of data at the 
7 TeV LHC. There is a combinatoric doubling of events, since 
incorrect leptonic W solutions are included. 



quired to reconstruct a W with \mjj — mw\ < 25 GeV. 
The effects of these cuts are shown in Table HI Even with- 
out the cuts, the signal is visible above the background 
for luminosity of order 1 fb~^, but this requires an ab- 
solute comparison to simulated backgrounds. The cuts 
allow backgrounds to be determined by a sideband anal- 
ysis, allowing for a robust discovery or exclusion. They 
are also required to reconstruct the h and the A. Note 
that these cuts do not depend on the h and A masses, and 
should be effective for a wide range of h and A masses. 

After these cuts, for an LHC luminosity of 1 fb^^ at 
7 TeV running we have approximately 10 signal events, 
corresponding to approximately 4.5(7 significance. With 
more data, one can reconstruct the h and A mass peaks. 
This is illustrated in Figs. [2] and [3J which are normalized 
to 1 fb~^, but include additional Monte Carlo statistics 
to show the shape of the peak. These plots include the 
incorrect reconstructions of the leptonic W. Additional 
techniques may be applied to reduce the effect of these 
incorrect solutions, but it is already clear that the mass 
can be reconstructed from these distributions. 

Conclusions — We have argued that models that have 



FIG. 3; Reconstructed A mass for 1 fb~^ of data at the 7 TeV 
LHC. There is a combinatoric doubling of events, incorrect 
leptonic W solutions are included. 



non-minimal Higgs sectors with heavier Higgs bosons are 
well- motivated, and generically lead to signals contain- 
ing 3 or more heavy standard model particles: f, Z, W, 
and b. We have advocated a simplified 2 Higgs doublet 
model to parameterize these signals. We have illustrated 
the possible searches with the process gg — >■ A^*'> -+ Zh 
followed either hy h^bb or WW/ZZ. For the WWZ fi- 
nal state we demonstrated that with as little as 1 fb^^ of 
data at the 7 TeV LHC this can lead to discovery of both 
the A and the h Higgs bosons. We hope that this work 
will motivate additional searches for final states involving 
3 or more heavy standard model particles. 
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